Two-color quantum memory in double A-media 



D. Viscor, 1 V. Ahufinger, 1 J. Mompart, 1 A. Zavatta, 2 ' 3 G. C. La Rocca, 4 and M. Artoni 2 ' 5 

1 Departament de Fisica, Universitat Autdnoma de Barcelona, E-08193 Bellaterra, Spain 
^European Laboratory for Nonlinear Spectroscopy (LENS), 
Via Nello Carrara 1, 1-50019 Sesto Fiorentino (Florence), Italy 
s ' Istituto Nazionale di Ottica, INO-CNR, L.go E. Fermi, 6, 1-50125, Florence, Italy 
Jf Scuola Normale Superiore and CNISM, Piazza dei Cavalieri, 56126 Pisa, Italy 
5 Department of Physics and Chemistry of Materials CNR-IDASC Sensor Lab, Brescia University, Brescia, Italy 

(Dated: January 16, 2013) 

We propose a quantum memory for a single-photon wave packet in a superposition of two different 
colors, i.e., two different frequency components, using the electromagnetically induced transparency 
technique in a double-A system. We examine a specific configuration in which the two frequency 
components are able to exchange energy through a four-wave mixing process as they propagate, 
so the state of the incident photon is recovered periodically at certain positions in the medium. 
We investigate the propagation dynamics as a function of the relative phase between the coupling 
beams and the input single-photon frequency components. Moreover, by considering time-dependent 
coupling beams, we numerically simulate the storage and retrieval of a two-frequency-component 
single-photon qubit. 



I. INTRODUCTION 

In recent years optical quantum memories had become 
the focus of an important research activity [IHZ] for be- 
ing one of the main ingredients for quantum informa- 
tion processing applications. In particular, in quantum 
information, the long distance transmission of photons 
(or flying qubits), which are the preferred information 
carriers, is limited by photon losses. Thus, transport- 
ing quantum states of light between different nodes of a 
quantum network requires the use of quantum repeaters 
[5J|!5], whose basic components are quantum memories. 
Therefore, quantum memories should be capable of stor- 
ing arbitrary quantum states of light for an arbitrarily 
long time and release them on demand and with high 
efficiency and fidelity pQ. 

Among the different methods for implementing a quan- 
tum memory, the approach based on electromagnetically 
induced transparency (EIT) [TDHT5] is one of the most 
used, allowing to store a single photon in solid state sys- 
tems for times > 1 s j!5j . This technique consists in slow- 
ing down a weak light pulse coupled to one transition of 
a A-type three-level system in the presence of a control 
field coupled to the other optical allowed transition. By 
adiabatically turning off the control field the light pulse 
is absorbed and mapped into the coherence between the 
ground states. Next, after a desired time which should be 
smaller than the decay time of the ground states coher- 
ence, the control field is turned on again and the initial 
light pulse is recovered. 

For the storage of a general photonic qubit, i.e., a sin- 
gle photon in an arbitrary superposition of two differ- 
ent components, more sophisticated schemes are needed 
[16H27| . In quantum communication with photons, the 
logical qubits can be encoded in several ways, for exam- 
ple, via polarization, time bin, path, phase, photon num- 
ber or even frequency encodings. In particular, several 
works have focused on the storage of photons with two 



frequency components using the EIT technique in reso- 
nant double-A media [23" 31! . Those proposals have been 
formulated mainly in the semiclassical regime. However, 
the storage of a two-color quantum entangled state would 
be interesting because it would have potential applica- 
tions in future quantum information networks, e.g., they 
could be used to link systems of different nature [32J [33] . 
One of the main issues regarding two-color memories, 
both in quantum and in classical approaches, is that the 
existence of a dark state in resonant double-A systems 
|34[ 135] together with the presence of four- wave mixing 
processes lead to a pulse matching of the frequency com- 
ponents [36] . This implies that the two input frequency 
modes can not be independently stored In particu- 
lar, only a specific combination of the two modes can be 
perfectly absorbed and recovered |3BJ [37] , whereas for 
an arbitrary two-frequency-mode input, part of the light 
will propagate transparently and part will be absorbed 
|24) . Nonetheless, it has been shown that the four- wave 
mixing processes arising in double-A media, which make 
difficult the implementation of a suitable quantum mem- 
ory, have interesting applications in frequency conversion 
of classical probe beams [3S1 135] j single-photon frequency 
conversion preserving the quantum coherence [35] , and in 
the possibility to combine or redistribute one or two pre- 
viously stored frequency modes [23] [27] [37], even with 
different relative intensities [21] • An interesting situation 
is found when one of the two A systems of the double-A 
configuration is far detuned from the one photon res- 
onance. In this case, considering the continuous wave 
regime, it has been shown that the total light intensity is 
weakly absorbed during the propagation [351 13H] , while 
the intensity of each mode oscillates sinusoidally with 
the optical length, being the energy transferred back and 
forth between the two probe beams. Later, a similar re- 
sult was obtained in the quantum regime [40] , where a 
single photon coupled initially to one of the transitions 
of the double-A system oscillates during propagation be- 
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tween the two frequency modes, thus creating a superpo- 
sition state at certain positions in the medium with high 
efficiency. 

In this work, we combine the usual EIT-based stor- 
age technique with the four-wave mixing properties of a 
double-A system to implement a quantum memory for 
single photons in an arbitrary superposition state of two 
frequencies. By solving the evolution equations of the 
single-photon frequency components we show that an ar- 
bitrary input superposition of two frequency modes can 
be recovered at certain positions of the medium, and that 
the relative phase between the coupling fields and the 
particular form of the input state play a crucial role in 
the propagation dynamics of the frequency components. 
Later, the storage and retrieval of the frequency super- 
position state is shown by the numerical integration of 
the system equations. 

The paper is organized as follows. In Sec. [Tl] we de- 
scribe the physical system and derive the equations that 
govern its evolution. In Sec. |III| we analytically solve 
the propagation equations of the incident single-photon 
frequency components and study several examples using 
different input superposition states and control fields pa- 
rameters. Next, in Sec. |IV[ numerical integration of the 
evolution equations of the system is performed to check 
the validity of the analytical approach, and the storage 
and retrieval of a particular input superposition state is 
presented. Finally, we summarize the results of this work 
and present the main conclusions in Sec. [V] 



II. THE MODEL 

We consider the physical system sketched in Fig. [T] 
where a single-photon wave packet in a superposition of 
two different frequency modes, of central frequencies ui^ 

and ujp 2 , and corresponding amplitudes E^y and E^ 2 , 
propagates through a system formed by A-type three- 
level atoms. Both frequency components interact with 
the left optical transition of the three-level atoms with a 
different detuning, 5 P \ or 5 P 2, being 5 P 2 far from the one 
photon resonance while 6 p i is close to resonance. We as- 
sume that the difference between the detunings is much 
larger than the spectral widths of the frequency compo- 
nents, such that there is no overlap between them. The 
other optical transition is driven by two strong coupling 
beams, of frequencies and tu^ 2 , tuned in two-photon 
resonance with the corresponding single-photon compo- 
nents, thus forming a double-A system. We consider that 
initially all the atoms are in the ground state 1 1) . The to- 
tal decay rate by spontaneous emission from the excited 
to the ground states is 72, and the decoherence rate of 
the ground states is denoted by 713. The total Hamil- 
tonian of the system is given by three contributions, the 
atomic {Ha), the field (Hp), and the interaction (Hj) 
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FIG. 1: (color online) . Double-A atomic scheme coupled with 
single-photon frequency components, E^ x and E^ 2 , and clas- 
sical field amplitudes, E c \ and E C 2, satisfying the two-photon 
resonance condition. The fields are detuned from the one pho- 
ton resonance with corresponding detunings 8 p i and 8 P 2, with 
S p i <C S P 2. All the atoms are initially in state 

Hamiltonians 

N 

H A = ^2J(W10-H + W 2 <722 + W3CT33), (1) 
i=l 

/oo 
SwpO^SwpdWp, ( 2 ) 
-OO 

N 

Hi = -E(^ 2 ^i ) ^ P + i + Mi2^ 2 1 ) -g+ 

i=l 

+ag ) hn° cl + ag ) hn c2 + H.c), (3) 

where the sums are over all the N atoms of the medium, 
the atomic population and coherence operators are of 
the form a vv = \v) (v\ and ai} = \v) (p\ , respec- 
tively, where v ^ p = {1,2,3} and j — {1,2} refers 
to the coherence generated by the mode uipj. The en- 
ergy of the atomic state \v) is given by hio u , h being the 
Planck's constant, /i 12 is the electric dipole moment of 
the |1) — 1 2) transition, and and a u are the cre- 
ation and annihilation field operators, respectively, for a 
frequency mode lj p . The Rabi frequencies of the clas- 
sical beams are denoted by Q° cj = ft ci e~^ ( *- z / c H^ 3 ; 
where Cl c j = I/X23 ] \E C j \ /fa, P23 is the dipole moment of 
the 1 3) — 1 2) transition and E c j the corresponding elec- 
tric field amplitude. The amplitudes of the quantum field 
operators read 

/oo 
eg^e+^'dwp, (4) 
-00 

(i) 

where c is the speed of light in vacuum and £tj p = 
e n [(lo p — uipj) / Awpj] , with e = \J~^^, £0 the electric 

permittivity in vacuum, OJ12 = ui\ — L02 the transition 
frequency between states |1) and |2), V the quantiza- 
tion volume, and n(w) a boxcar function of width Au> p j, 
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centered at w pj -. We assume Auj p j much larger than the 
spectral width of the corresponding frequency component 



Suipj, but not enough to overlap with the other one, i.e. 

| W p2 ~~ "'pi | 



w2i I > Aw. 



To find the evolution equations of the single-photon 
frequency components we adopt the procedure and for- 
malism from Ref. [TTJ. The initial state of the system 
has the form 



fofc(t->-oo)) - J dw pl /W(-oo)at pi |o,0) p |l) 

+ / dup2f™{-ao)3i J0,0) P |1). (5) 



where u> p j G [w. 



PJ 



■Au pj /2,u)° 



+ Aoj p j/2]. Here we have 
used the notation |ni, 712)3, where ni and n 2 are the 
number of photons in modes w pl and uj p2 , respectively, 
and v denotes the atomic state. Tracing out the atomic 
part, the first and second terms in the right-hand side of 
Eq. ([5| correspond to the initial state of each frequency 

component of the single photon, and fujj— 00) are the 
envelope functions of the wave packet, which have a nar- 
row peak at 0;° • . We assume that they are spectrally sep- 
arated enough such that their overlap is negligible, and 



/ffi(-oo) 



■J dui : 



P2 



/$(-«>) 



must satisfy J dui p i 

1. So, the photon is initially in an arbitrary superpo- 
sition of the two frequency components, and the atoms 
are all in the ground state |1). Next, we assume that the 
general form of the state of the system at any time is 



|V(*)) = l^l(*)) + l^(*)> + l^3(t)), 



(6) 



where the first, the second and the third terms corre- 
spond to the excitation being either in one of the pho- 
tonic modes, in the atomic state |2), or in state |3), re- 
spectively. Their explicit forms are 



the Schrodinger equation and apply (1,0| P (1|, (0, l| p (l|, 
(0,0|p(2|, and (0,0| p (3|, obtaining: 



u 'pj fjpj 


(t) 




~ gm% 


M12 


f *» V if!& 


U3g(t) - 





Nebj(t)e- 1 - 



(10) 
(11) 



(12) 

and integrating 



Next, multiplying Eq. (10 1 by e luJpjZ / c 
over oj v j we obtain the propagation equation for the 
quantum field amplitudes 

1 



d t + d z ) £j(z,t) = iK 12 Pj(z,t), 



(13) 



where we have defined £j(z,t)e 



and k X2 = 



ILUpj z/c 
2 2 



Joj° At-z/c) 



and (12) read 



&(z,f) = bj{t)e* 
With these definitions, Eqs. |IT| 



d t g(z,t) 



iA p j0j(z,t) +i£j(z,t) + ig(z,t)Cl c j, (14) 
T%lPi(z, t) + Sf c2 i3 2 (z, i)] - 7l3 ff(z, t), 

(15) 



where we have added phenomenologically the ground- 
states decoherence 713 in Eq. ( 15 ) and the spontaneous 



emission from the excited level 72 in Eq. ( 14 ) through the 

complex detuning A p j — S P j— 272/2, being S P j = oj^—uj 2 . 

Moreover, we have defined fl C j = 

sumed degenerate ground states 

w°. - 
pj 

huj\ = 0. 



n cj -e i we have as- 
i.e., u>i 

i^cji an( l we nave chosen the energy origin at 



W3 and 



|V>i(*)) - J ^ P xf^{t)ai pl %Q) p \l) 

+ J du p2 f^{t)aijQ,Q) p \l), (7) 

N 

llM*)> = ^(6i(i)^ 1 i ) |0,0) p |l)+6 2 (t),7( 2 1 ) |0,0)p|l)), 

(8) 

N 

llM*)> = ^ff(*)^3i|0,0>p|l), (9) 



where bi{t) and b 2 (t) are the probability amplitudes of 
exciting one atom to the state |2) through modes ui pl 
and ojp 2 , respectively, and g(t) is the probability ampli- 
tude of transferring the population to state |3) via a two- 
photon process. Those functions together with f!^ p { (t) 
(2) 

and /w p2 (i), give a complete description of the state of 
the system. In order to find their evolution, we insert 
the general form of the state of the system [Eq. ([6|] into 



III. SOLUTIONS OF THE EVOLUTION 
EQUATIONS 

The equations describing the propagation of the single- 
photon frequency components in the double-A system 
can be solved by using the adiabatic approximation for 
Eq. (14 1, i.e., dt/3j(z,t) ~ 0, and changing from temporal 
to frequency domain by applying the Fourier transform to 
our system equations. Next, inserting Eqs. (14 1 and (15) 



into the Fourier transformed Eq. (13), a linear system of 
partial differential equations for the quantum field am- 
plitudes in the frequency domain, £j(z,u>), is obtained. 
This can be solved, leading to 



£, 



[z,u) = ^(O,^^ 2 ( e*"/*. + J_y! e W» e *» 



in 



n 



(16) 
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where j,l = {1,2} and j ^ I, £j(Q,u)) is the boundary 
condition for the spectral envelope of the frequency 

n lo ' 2 
ponent centered at w^-, a 



com- 



1 

v a 
1 

"6 



«12 



-^i2 J -g- and 



Kl2 |^cl| 2 |^ c2 | 2 (Api - A p2 ) 2 



Ifil 



D 2 



(17) 



(18) 



with D = A pl \fl 

\n 



c2 



c2l "I" A p2 |^ c i| z and |f2| z = |fi c i| z + 
In Eq. (|16), we have assumed that the decoher- 



ence time of the ground states is much larger than the 
time needed to store and retrieve the single photon, thus 
713 ~ 0. Moreover, we have approximated the expo- 
nents as linear functions of w by Taylor expansion up 
to first order, and we have considered the coefficients 
independent of u>, as done in Ref. [ID]. With these ap- 
proximations, and considering a small decay from the 



excited level 72 <§C (<5 p i |n c2 | + 8 p2 



the 



verse Fourier transform of the field can be performed an- 



alytically. Then, it can be seen that in general each of 
the components of the frequency superposition will split 
in two different parts, each one propagating with a dif- 
ferent velocity given by Eqs. (17) and (18). Assuming 
|fici| — |0 C 2|, dpi — 0, and 5 P 2 3> 72, the velocities 
for the frequency components are approximately equal 
Vb ~ v a = v, and hence the inverse Fourier transform of 



Eq. ( 16 ) can be rewritten as 



£j(z,t) = 



(*-;)< 



1 + e 1 



+£1 [0,t- -J e 1 ^' (l-e iQ2 )J , (19) 



where (fiji = 4>j — <j>i is the phase difference between the 
coupling fields. Note that a now reduces to a ~ + 

i«i 2 f^ 02]. 

p2 

To obtain the intensity of each component of 
the single-photon frequency superposition we calculate 
£j(z,t)£j(z,t) using Eq. (19): 



\£ 3 (z,t)\ 2 = 



+2Re 



^ (* ~ ~) I { 1 + e_2Im(Q) " + 2 cos [Re(a)z] e - Im(Q)z } + 

° (* ~ ") 2 { 1 + e ~ 2lm(a)z - 2cos [Re(a)z] e- Im < Q ^} + 

e mi e i<l>ii |i _ e -2im(a)* +2 isin[Re(a)z] e - Im ( Q ) 2 |) 



5° it- 



(20) 



where y>,-/ = (pj — ipi is the phase difference between the single-photon frequency components at the input z — 0. In 
Eq. (20) the time evolution appears only in the boundary conditions £® (t — ^) = £j(0,t — ^). This means that the 
single-photon wave packet keeps its shape, but it is drifted in time a quantity t c — z/v, which depends on the velocity 



v defined in Eq. (17). Moreover, it can be seen that, while the single photon propagates, the intensities of the two 



frequency components exhibit complementary oscillations with a rate that depends on Re(a). Note that the decaying 
terms Im(a) in Eq. (20) are due to spontaneous emission from the excited level. An interesting case is found when 
one considers a symmetric superposition state at the input, i.e., |£? (t)\ = (t)\. In this situation, Eq. (20) takes 
the form 



\£ 3 ( Zl t)\ 2 = \^(t- Z -) 



2 / l _l_ e -2Im(a)z 



'2 g— 2Im(a)j2 

+ <j cos (<j} jt + tpji) ( ) - sin {4>ji + <p jt ) sin [Re(a)z\ e^ Im(a)z 



(21) 



In this case, when <pji + tpji = 0, there is no oscillation between the frequency components during the propagation 

and the intensity of each single-photon frequency component is perfectly transmitted, i.e., \£j(z,t)\ 2 = \£® (t— f)| 2 - 
This can be interpreted by considering that, through a four-wave mixing process mediated by the coupling beams, 
the energy going from the first to the second component is compensated by the energy transfer from the second 
component to the first one. 
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Analogously, we find the relative phase between the two frequency components using 

/ / Z\ / Z\ ( J I „-2Im(u)z 

£ 3 (z,t)£^(z,t) = iS°(t--) £?[t--) | cos foi - +isin(^ 7 -^)e _Im(Q)2 cos[Re(a)^] 



\£° (t 

+i sin [Re(a)z] e ~ Im ( Q ) 2 



1 _ e -2 M a), |gO^„ f )|V|gO( t _ g )| 

2 2 



(22) 



From this expression, we observe that in general the 
phase between the two components, arg [£j(z,t)£f(z,t)], 
will oscillate in a more involved way than the intensity, 
[Eq. (20 1]. In particular, we observe that only when the 
imaginary part of the outermost parentheses in Eq. ( 22 ) 
vanishes the phase will be independent of z. 

In what follows, by evaluating the analytical expres- 
sions obtained in Eqs. (20 1 and (22 1, we discuss different 



propagation examples of the two frequency components; 
see Figs. [2] and [3j We change to a reference frame fixed 
at the peak of the single-photon pulse (t c — z/v), so we 
need only to show the variation on the spatial dimension 
z. In Figs.[2ja) and|3^a) we plot the normalized intensity 
of each of the two frequency components, 



I 3 {z) 



i£?(Mc)r 



\£Ut c )\ 2 + \£Utc)\ 2 ' 



(23) 



whereas in Figs.J^b) and^b), the relative phase between 
them, 



$jl{z) = arg [Sj(z,t c )£i(z,t e )], 



(24) 



is shown. In all the figures the different line styles cor- 
respond to different sets of parameters (see the caption) , 
while the black and gray lines both in Figs.[2ja) and[3^a) 
correspond to the intensity of the frequency components 
and ujp 2 , respectively. We have taken Im(a) = 0, a 
fact that is well justified from the assumption 5 P 2 3> 72 
made in Eq. (19). These figures are useful to show that 



the behavior of the two components during the propa- 
gation depends completely on the specific state at the 
entrance of the medium and the phase difference of the 
coupling beams. For instance, the different line styles 
in Fig. [2] correspond to different initial intensities of the 
frequency components, while the relative phases between 
the frequency components and coupling beams are fixed. 
We observe that the different initial superposition states 
lead to intensity oscillations with different amplitudes 
and shifted by different amounts. Further examples are 
shown in Fig. [3j where the input intensities are equal for 
the two frequency modes, and the relative phases between 
them and between the coupling beams are changed. We 
observe in Fig. [3]ja) that opposite behaviors for the in- 
tensity of a given mode are obtained just by properly 
changing the relative phase of the coupling beams (solid 



and dashed lines). Moreover, note that the case shown 
with dotted lines, i.e., I p x(0) — I P 2(0) and cf>i2 = fi2 = 0, 



corresponds to the situation discussed after Eq. (21), in 



which the photon state does not evolve during propaga- 
tion. 

As a general conclusion from Figs. [2] and |3j we observe 
that the more different the intensities of the frequency 
components, the largest the variation in their relative 
phase, and vice versa. We also observe that the relative 
phase oscillates around the value cj>ij. Moreover, the most 
remarkable fact is that the frequency of the oscillation, 



both in the intensity [Eq. ( 23 1] and in the phase [Eq. ( 24 )] 
is determined only by Re(a) ~ — 2ki2/<5 P 2- This means 
that by properly choosing the coupling parameter K12 
and the detuning 5 P 2, one can recover at the output of the 
medium, z = L, the initially injected state with an ideally 
perfect fidelity, i.e., £j(L,t) — £j(t — t c ) for Re(a)L = 
27rn, with n € Z. 



IV. NUMERICAL ANALYSIS 

In this section we demonstrate the validity of the ap- 
proximations made in the analytical approach by numer- 
ically integrating Eqs. ( 13 )-( 15 1. Moreover, we show the 



possibility of storing and retrieving a single photon in an 
arbitrary superposition state of two frequency compo- 
nents using time-dependent coupling fields. To simulate 
the pulse propagation in time and space, a bidimensional 
grid for each variable is created with a spacing in the 
z dimension small enough to ensure the convergence of 
the results. The steps for the numerical protocol are the 
following: First, the temporal evolution of the medium 
variables is obtained from the incident (z = 0) field com- 
ponents, which are assumed to have Gaussian profiles of 
temporal width t = 25 ns and centered at t c — 3.5r, 
using a Runge-Kutta integrating method. Next, the field 
at the adjacent spatial point is determined with a finite 
difference method, using the preceding obtained values. 
Finally the previous steps are repeatedly performed un- 
til the whole grid is filled. For the medium we take a 
length of L = 0.1 m, 72 r = 0.16, 7 i 3 r = 1.6 x 10~ 5 , 
and k\ 2 tL ~ 500, while the detunings are S p \t = and 
S p2 r = 160. 

On the one hand, an example of the propagation of 
the two single-photon frequency components is shown in 
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FIG. 2: (a) Normalized intensities of the single-photon fre- 
quency components 7i (black) and I2 (gray), and (b) the rel- 
ative phase between them for Re(a) = 2n/L, with L — 0.1 
(a.u.) being the length of the medium, and Im(a) = 0. 
Solid lines: 7i(0) = 0.99, 0i2 = 0, 1^12 — n/4; dashed lines: 
7i(0) = 0.85, 0i2 = 0, ip 12 = tt/4; dotted lines: 7i(0) = 0.70, 
0i2 = 0, tpia = 7r/4. 



FIG. 3: (a) Normalized intensities of the single-photon fre- 
quency components I\ (black) and I2 (gray), and (b) the rel- 
ative phase between them for Re(a) = 2ir/L, with L — 0.1 
(a.u.) being the length of the medium, and Im(a) = 0. 



Solid lines: 7i(0) = 0.5, 012 

h(0) = 0.5, 012 = 7T/2, <fi 12 = 
012 = 0, (fl2 = 0. 



= 0, (fi2 = tt/3; dashed lines: 
tt/3; dotted lines: 7i(0) = 0.5, 



Fig. [4j where the normalized intensity of both compo- 
nents, I\ (a) and I2 (b), is shown as a function of po- 
sition and time for constant coupling Rabi frequencies 
I O CJ - 1 t = |fi c ; t = 18, and a phase between them of 
012 = 0. The peak amplitudes of the frequency com- 
ponents for the injected single photon are |£°(t c )| r = 

1.3 x 1(T 3 and \£${t c )\ = ^0.3/0.7 |£?(t c ) | , with a rel- 
ative phase (pi2 = 7r/4 between them. Note that those 
parameters correspond to the case represented in Fig. [2] 
with dotted lines. As we observe, the intensities for the 
two frequency components exhibit complementary oscil- 
lations with a spatial period of ~ 0.125L. Moreover, the 
displacement of the peak allows to estimate a propaga- 
tion velocity of ~ 10 6 m/s. Using the model derived in 
the previous section, the values obtained for the oscilla- 
tion period and the velocity are ~ 0.1L and ~ 4.5 x 10 6 
m/s, respectively. Thus, the numerical simulations are 
in agreement with the analytical results. The phase be- 
tween the components at the peak of the pulse, $ji(z), is 
plotted in Fig. |3]Jc) as a function of z. We observe that 
the behavior of the phase is in good agreement with the 
analytical result [see dotted line in Fig. |2fb)]. 

On the other hand, Fig. [5] shows a particular exam- 
ple of the storage and retrieval process, using temporal 
profiles of the coupling beams of the form 

U' jt) = {2 - tanh [a(t - ti)] + tanh [a(t-t 2 )}}, 
2 




0.0 0.1 0.2 0.3 0.4 0.5 



z/L 

FIG. 4: (color online). Normalized intensities of the single 
pulse frequency components, 7i (a) and 72 (b), as a function 
of normalized position and time, and (c) the phase between 
the frequency components at the pulse peak as a function of 
normalized position. The parameters correspond to the case 
represented with dotted lines in Fig. [2] 
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V = T) Abs??Rct ! with the absorption ryAbs and retrieval r)R Ct 
efficiencies being 



FIG. 5: (color online). Temporal profile of the coupling beams 
(a) and normalized intensities of the single pulse frequency 
components, I\ (b) and I2 (c), as a function of normalized 
position and time. The parameters correspond to the case 
represented with dotted lines in Fig. [3] 



with |r2 cl |T = |fi c2 |r = 18, a = 0.5r _1 , ti = 2t c 
and ti — 6t c [see Fig. [5^a)], and a phase difference be- 
tween the coupling fields of $12 = 0. In Figs. [5jb) and 
(c) the normalized intensity of frequency components 
1 1 and I 2 , respectively, is shown as a function of po- 
sition and time. In the example we have taken equal 
amplitudes for the two components of the input state, 
|^i {tc) I T = I f 2 (*c ) I t = 1.3 x 1CP 3 , and an initial phase 
difference ip%2 — between them, in such a way that 
the chosen parameters correspond to the situation with 
constant coupling fields represented with dotted lines in 
Fig. [3] Figure [5] shows an example of how the superposi- 
tion state can be stored and recovered by appropriately 
varying in time the coupling fields. Here, the storage 
time corresponds approximately to ^2 — ^i = 0.35 /zs, and 
it could be extended in principle to times of the order 
°f I/713 (~ 1-5 ms for the parameters considered). The 
behavior of the intensities for each component coincides 
with the predictions of the theoretical model. We have 
checked that the total pulse area is almost conserved al- 
though the pulse spreads during propagation. The phase 
between the frequency components (not shown in the fig- 
ure) keeps an approximately constant value of <&ji{z) — 
during the whole storage and retrieval process, as ex- 
pected from the dotted line in Fig. [3^b) . To characterize 
the memory performance, the efficiency of the storage 
and the retrieval processes, and the fidelity of the re- 
covered superposition state have been calculated. On 
the one hand, we define the performance efficiency as 



?7Abs 



f t f (|£ 1 (L,t)| a + |£2(£,t)| 2 )dt 
1 )—. ^ — , (25) 

ill 1 (ism 
ft! /2 (W(t)\ 2 + \m)t 



(26) 



where the interval t = 0, tf — 30t is the integration 
time. Computing these expressions with the data ob- 
tained in the simulation shown in Fig. [5j the absorp- 
tion and retrieval efficiencies are ??Abs = 99.78% and 
^Rot = 91.21%, respectively. Thus, the total efficiency 
is i) — 91.01%. On the other hand, the conditional fi- 
delity is defined as F c — \(ipi n \ Vw)| 2 , where we take as 
input and output states 



St ,/a \%(!)\'*t 



/ t * //2 (i^(*)i 2 +i^ (*)i a )* 

e<^)fl f/2 \£$(t)\ 2 dt 



IVWt) 



jy /2 (is^tf+ML^dt 

^St^ ML,t)\ 2 dt 
/// /a (\£,{L,ttf + \£2(L,t)\ 2 )dt 



|1,0)J1) + 



|0,1)J1>, (27) 



|l,0) p |l) + 



|0,1>J1) 



(28) 



respectively, with (ipy^) = f t * arg [£i(Q, (0, t)]dt and 
(tplf) = j/^axg^Z, Therefore, the cal- 

culated conditional fidelity for the case shown in Fig. [5] 
is F c = 99.69%. 



V. CONCLUSIONS 

In this work we have studied the propagation of a sin- 
gle photon, in an arbitrary superposition of two different 
frequency components, through a double- A medium. In 
this particular configuration the intensities of the two fre- 
quency components exhibit complementary periodic os- 
cillations as they propagate. These propagation effects 
have been used in combination with the light storage 
technique based on EIT to implement a quantum mem- 
ory for frequency encoded single-photon qubits. We have 
studied analytically the dependence of the relative phase 
between the coupling fields and the input qubit state 
in the propagation dynamics. Moreover we have shown 
that, at certain positions in the medium, the initial single 
photon, which can be in any desired frequency superposi- 
tion state of the two frequency components, is recovered. 
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The numerical results, obtained by numerically integrat- 
ing the evolution equations of the system, are in good 
agreement with the analytical solutions and thus the va- 
lidity of the analytical approach has been confirmed. Fi- 
nally, the storage and retrieval of a single-photon state in 
an arbitrary superposition of two frequency components 
has been shown numerically by turning off and on the 
coupling fields during the propagation of the single pho- 
ton. For the specific choice of parameters we adopt here, 
the results demonstrate an efficient quantum memory for 
high-fidelity storage and retrieval of a frequency encoded 
single-photon qubit. 
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